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Abstract 

An  integral  transform  analysis  of  the  static  and  time-harmonic  scattering  of  the  scalar 
(two-dimensional)  field  radiated  by  a  line  source  in  the  vicinity  of  a  penetrable  and 
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density  contrast  wedge  also  gives  analytical  solutions  of  the  Helmholtz  equation,  in  a  form 
suitable  for  the  asymptotic  extraction  of  the  high-frequency  ray  components.  Similarly, 
the  transform  analysis  of  the  impedance  boundary  wedge  results  in  a  difference  equation 
to  be  solved  in  the  transform  variable.  A  special  inhomogeneous  surface  impedance  yields 
purely  algebraic  equations  for  the  transforms,  which  can  be  solved  in  closed-form  and  in¬ 
verted.  At  each  stage,  the  mathematics  and  discussion  are  guided  by  the  physics  of  both 
the  acoustic  and  equivalent  electromagnetic  scattering  problems. 
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I.  Introduction 


The  intentional  “or”  in  the  title  of  this  report  emphasizes  the  exact,  mathematical  equiv¬ 
alence  between  the  acoustic  and  electromagnetic  scattering  problems  for  a  two-dimensional, 
penetrable  wedge.  A  z— directed,  time-harmonic  line  source  at  the  transverse  (x,  y)  posi¬ 
tion  f '  is  external  and  parallel  to  the  edge  of  an  infinite  wedge  of  included  angle  2a.  Both 
the  exterior  (region  1)  and  interior  (region  2)  are  composed  of  simple  media  having  linear, 
isotropic,  and  homogeneous  properties.  See  Fig.  1  of  page  7,  for  example. 

In  the  acoustic  problem  for  the  velocity /density  contrast  wedge,  which  is  generally 
called  the  transmission  problem  in  the  mathematics  literature,  the  media  ( j  =  1  or  2) 
constitutive  parameters  are  conveniently  taken  to  be  the  ambient  density  pQi  and  the  wave 
speed  Cj.  In  the  electromagnetic  problem,  each  medium  is  characterized  by  its  electric 
permittivity  Cj  and  magnetic  permeability  pj ,  which  naturally  appear  together  as  the  wave 
speed  Cj  =  l/y/pjej  and  intrinsic  impedance  rjj  =  y/pjjfj-  The  electromagnetic  scatterer 
is  referred  to  as  the  dielectric  wedge,  even  though  p  as  well  as  e  can  independently  differ 
in  regions  1  and  2. 

At  the  constant  angular  frequency  u>,  associated  with  each  unbounded  medium  is  the 
wavenumber  kj  =  w/cj.  The  reduced  wave  or  Helmholtz  equation 

(VJ  +  *?)  V-(r-)  =  -«(?-  ?') 


governs  the  total  field  in  each  of  these  three  different  physical  situations: 

1.  Acoustic  line  source,  with  ip(r)  the  scalar  pressure  field.  Boundary  conditions  at 
the  material  interfaces  are  continuity  of  the  pressure  and  normal  velocity, 


ip\  =  and 


1  1  d  ,  115, 

- *7^1  = - 

Poi  r  0<p  po2  r  o<f> 


2.  Magnetic  line  source,  where  the  scalar  field  of  interest  is  the  single  z— directed 
component  of  magnetic  field  Hz{f)  =  ip(f)  in  this  often-called  case  of  TE  or 
transverse-electric  polarization.  A  unique  solution  to  Maxwell’s  equations  requires 
continuity  of  tangential  magnetic  and  electric  fields, 


ip  i  =  V*2  and 


1 1  a  .  i  i  a  , 

- njV’i  = - ^7^2 

6i  r  0<p  €2  r  0<p 


3.  Electric  line  source,  where  the  scalar  field  of  interest  is  the  single  z—  directed  com¬ 
ponent  of  electric  field  Ex(r)  =  ip{r)  in  this  often-called  case  of  TM  or  transverse- 
magnetic  polarization.  Boundary  conditions  of  continuity  of  tangential  electric  and 
magnetic  fields  are 


=  i>2  and 


1 1  a  ,  i  i  a  . 

- = - Tl'P  2 

pi  r  0<p  pzr  o<t> 
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All  three  of  these  different  physical  problems  are  governed  by  the  same  differential  equation, 
with  different  physical  origins  for  the  constants  appearing  in  the  mathematically  identical 
boundary  conditions.  Therefore,  this  acoustic  and  electromagnetic  problem  is  properly 
treated  as  one  boundary  value  problem,  and  liberal  use  is  made  of  the  archival  contributions 
of  acousticians,  electricians,  and  wave  mechanics. 

The  perfectly  soft  (Dirichlet  or  TM  boundary  condition)  or  perfectly  hard  (Neumann 
or  TE  boundary  condition)  impenetrable  wedges  are  historically  solved  by  separation  of 
variables.  Familiarity  with  these  classical  solutions  is  important  for  guidance  and  sanity 
checks  of  any  approximate  methods  applied  to  the  penetrable  wedge.  The  Kontorovich- 
Lebedev  transform  (a  spectral  superposition  of  cylindrical  waves)  of  D.  S.  Jones  [49]  is 
extended  in  Sections  III  and  IV  to  handle  the  density  contrast  and  impedance  boundary 
wedges,  respectively.  The  analysis  of  the  impedance  boundary  wedge  results  in  a  difference 
equation  in  the  transform  variable,  which  is  interesting  in  view  of  the  apparently  unrelated 
difference  equation  that  appears  in  a  1959  paper  by  W.  E.  Williams  [7]. 

The  analytic  solution  of  the  Leontovich  impedance  boundary  conditions  by  Maliuzhinets 
[6]  in  1958  is  still  being  asymptotically  evaluated  [9-12,25-26]  to  derive  high-frequency 
diffraction  coefficients  in  the  spirit  of  the  geometric  theory  of  diffraction  (GTD)  or  related 
wedge  assemblage  [41].  Heuristic  diffraction  coefficients  [27]  consisting  of  the  product 
of  the  accepted  forms  for  perfect  conductors  and  plane-wave  reflection  coefficients  at  a 
dielectric  interface  are  absolutely  invalid.  A  seemingly  ad-hoc  correction  to  the  physical 
optics  or  Kirchhoff  approximation  [24]  gives  results  that  agree  with  numerical  experiments 
in  some  instances,  but  offer  little  hope  of  getting  closer  to  the  important  wave  physics. 
Obviously  then,  the  prospect  of  rigorously  deriving  the  high-frequency  ray  behavior  is  a 
prime  motivation  for  the  present  research  effort. 

Although  it  may  initially  appear  strange  to  the  uninitiated,  the  static  solution  is  a 
critical  component  of  any  dynamic  scattering  problem,  especially  with  regard  to  edge 
conditions  and  behavior  in  source  regions.  Section  II  is  a  complete  and  exact  solution  to 
the  static  dielectric  wedge  problem.  Based  on  literature  searches  and  discussion  with  E. 
Marx  (author  of  [20-22,  28]),  it  is  presently  believed  that  this  static  solution  is  new.  The 
anticipated  peer  review  of  a  refereed  journal  submission  based  on  Section  II  should  confirm 
or  refute  this.  It  should  be  noted  that  the  special  case  with  the  line  source  lying  on  the 
x— axis  (the  wedge  bisector)  does  appear  in  the  Russian  text  [53].  Furthermore,  Smythe 
[50]  presents  a  formal  treatment  of  the  static  dielectric  wedge,  which  is  unfortunately  flawed 
with  divergent  integral  representations.  In  any  event,  the  exact  static  solution  of  Section 
II  is  now  available  to  confidently  answer  any  questions  about  the  near-in  edge  behavior. 
The  static  problem  is  solved  using  the  Mellin  transform,  which  is  itself  the  static  limit  of 
the  Kontorovich-Lebedev  transform  employed  in  Sections  III  and  IV. 

It  is  important  to  address  the  literature  on  the  penetrable  wedge,  and  to  put  this 
research  in  historical  perspective.  Rawlins  [1]  constructs  a  perturbation  solution  of  a 
volume  Fredholm  integral  equation  of  the  first  kind  (actually  a  surface  equation  in  R2), 
which  converges  for  low  values  of  the  relative  permittivity  or  dielectric  constant  er  =  e2/ci . 
Several  authors,  including  Kleinman  and  Martin  [4],  prove  convergence  for  1  <  <  2, 

and  divergence  of  the  Neumann  iteration  series  for  y/e7  >  2.  Except  for  largely  numerical 
schemes,  the  1979  observation  “ . . .  little  progress  has  been  made  with  solving  the  canonical 
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problem  ...  ”  by  D.  S.  Jones  [48]  has  remained  true.  Fortunately,  the  progress  and  solid 
foundation  of  the  present  research  is  encouraging,  especially  with  regard  to  the  novel 
solution  of  the  impedance  boundary  wedge. 

In  competition  with  the  integral  transforms  that  have  become  the  method  of  choice  in 
this  first  phase  of  the  project,  surface  integral  equations  are  applicable  to  the  transmission 
problem.  Colton  and  Kress  [50]  give  one  of  the  most  complete,  and  therefore  theoretical, 
accounts  of  coupled  Fredholm  integral  equations  for  equivalent  surface  distributions.  The 
infinite  extent  of  the  surface  may  appear  a  little  discouraging,  but  perhaps  no  more  so 
than  the  infinite  domains  in  the  integral  transforms.  Davey  [5]  gets  some  mileage  out  of 
a  physically  intuitive  subtraction  of  the  expected  far  (r  — »  oo)  surface  behavior,  which 
is  to  first  order  the  behavior  at  an  infinite,  planar  boundary.  Glisson  [19],  Marx  [20], 
and  Kleinman  and  Martin  [4]  combine  the  pair  of  coupled  integral  equations  into  a  double 
integral  operator  of  a  single  unknown  surface  function,  which  they  all  claim  is  preferable  to 
two  unknown  functions.  Two  coupled,  single  integral  operators  could  in  fact  be  analytically 
and  computationally  preferable,  especially  in  view  of  the  “hypersingular”  kernel  in  the 
double  operator. 

The  substantial  contribution  by  Chu  [45],  in  1989,  for  the  impulse  response  of  the  den¬ 
sity  contrast  wedge  can  serve  as  a  future  check  on  the  time-harmonic  results  of  Section  III, 
after  a  suitable  Fourier  transformation  between  time  and  frequency.  A  three-dimensional 
plane  wave  reflection  coefficient  using  a  finite  number  of  images  is  assembled  by  Deane  and 
Tindle  [46],  in  one  of  many  JASA  papers  that  underscore  the  applicability  of  propagation 
in  wedge-shaped  regions  to  shallow  ocean  acoustics.  Westwood  [47]  employs  an  approxi¬ 
mate  summation  (integral)  of  complex  ray  contributions  in  the  spirit  of  Felsen,  and  then 
computes  some  wideband  transient  responses  via  the  Fourier  transform. 

The  research  summarized  in  this  report  is  a  first-principles  effort  toward  a  fundamental 
understanding  of  the  static  and  time-harmonic  excitation  of  the  dielectric  (penetrable) 
wedge.  Arbitrary  transient  waveforms  are  not  addressed  in  this  phase  of  the  research,  as 
a  Fourier  transform  of  the  time-harmonic  response  is  the  traditional  and  logical  approach 
to  time  domain  scattering.  The  case  of  plane  wave  excitation  is  obtained  by  letting  the 
source  recede  to  infinity,  i.e.  r'  — »  oo.  Each  section  is  mostly  self-contained,  although 
frequent  reference  is  made  to  the  details  of  the  completed  static  solution  of  Section  II.  A 
concise  summary  and  suggestions  for  the  most  fruitful  avenues  to  pursue  for  the  coming 
year  appear  in  Section  V. 
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II.  Exact  Static  Solution  for  the  Dielectric  Wedge 


Fig.  1.  Dielectric  Wedge  and  Static  Line  Charge 

The  static  ( Jb  =  0)  excitation  of  the  penetrable  wedge  is  treated  in  the  electrical  parlance, 
where  the  line  charge  of  unit  lineal  density  is  located  at  the  source  coordinates  (r',  <f>')  of 
Fig.  1.  The  permittivity  of  the  wedge  of  angle  2a  is  ej,  which  is  surrounded  by  a  medium 
with  permittivity  ei .  Consistent  with  the  other  chapters  of  this  report,  all  geometry  is  of 
infinite  extent  and  invariant  in  the  z— dimension,  which  restricts  the  physical  domain  to 
R3. 

The  irrotational  electrostatic  field  is  uniquely  characterized  by  the  scalar  potential 
which  is  a  solution  of  Poisson’s  equation  subject  to  appropriate  boundary  conditions. 
Denote  by  4>\(r)  the  potential  field  in  the  external  region  where  the  source  is 

VViM)  =  6(r  -  r')6(<fi  -  <f>')  (a  <  <f>  <  2*  -  a)  (1) 


and  let  ipi(r)  be  the  source-free  field  inside  the  wedge 

V2^(r,  <t>)  ~  0  <  <*)•  (2) 

Boundary  conditions  at  the  material  interfaces  <f>  =  ±a  for  this  scalar  potential  are  con¬ 
tinuity  of  ip  (from  continuity  of  the  tangential  electric  field)  and  continuity  of  the  normal 
electric  flux  density  drpj  / dn  (absence  of  any  free  surface  charge). 

In  order  to  simplify  the  ensuing  analysis,  it  is  expedient  to  decompose  the  desired 
solution  rp(x,  y)  for  the  boundary  value  problem  of  Fig.  1  into  its  odd  rp°(xty)  and  even 
ipe(x  ,y)  components 


rp{x,  ±y)  =  i  tye(*,  y)  ±  rp°(x,  y)]  (y  >  0)  (3) 

with  respect  to  the  x— axis  which  bisects  the  wedge.  Since  ip°(x,  y)  is  an  odd  function 
of  y,  it  vanishes  on  the  y  =  0  plane  and  is  therefore  the  solution  in  the  upper  half-space 
y  >  0  for  the  problem  having  a  soft  bisecting  plane  (Dirichlet  boundary  condition).  This  is 
equivalent  to  an  out-of-phase  image  source.  Similarly,  rpe(x,  y)  is  the  solution  for  the  hard 
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Fig.  2.  Odd  and  Even  Symmetry  Components 


bisecting  plane  (Neumann  boundary  condition),  which  sustains  an  in-phase  image  source. 
This  symmetry  is  depicted  in  Fig.  2. 

A.  Mellin  Transform  for  the  Case  of  Odd  Symmetry. 

The  source  coordinate  ^  =  <j>'  divides  region  (1)  into  two  source-free  regions,  resulting 
in  a  total  of  three  subregions  to  consider  for  the  half-space  above  the  perfectly  soft  plane: 


V’aO'j  <t>),  0  <<f><a 

y)  =  '  V»f  (r,  4>),  a  (4) 

,  <f>'  <  $  <  it. 

The  Dirichlet  boundary  condition  on  the  soft  plane  is 

(r,  0)  =  0  (5) 

ft  (»*,  tt)  =  0  (6) 

and  at  the  material  interface 

a)  =  V’f  (»*,  a)  (7) 

e2^'2(r’°)  ~  ei0?^r,a)'  ^ 

Consistent  with  the  conventional  Green’s  function  ansatz,  the  potential  is  continuous 
everywhere  across  the  source  plane  <f>  =  ft 

(r>  ^')  =  V'i (9) 

while  the  discontinuity  in  normal  derivative 

jyK  <r-  *">  -  ^  (*•  *')  =  77*<r  -  r'>  (10) 


results  from  integrating  (1)  from  4>'~  to  4>'+- 

The  Mellin  transform  [52]  is  used  to  transform  the  radial  variable  r  to  a  complex  variable 
s,  whereupon  the  remaining  differential  equation  with  boundary  conditions  in  <f>  is  solved  in 


closed-form.  The  physical  solution  tp(r,  <f>)  is  recovered  by  careful  evaluation  of  the  inverse 
Mellin  transform.  The  Mellin  transform  of  /(r)  is 

OO 

F(s)  =  M{f(r);i)  =  J  T'~'f(r)dr  (11) 

0 

and  has 

e+ioo 

/(r)  =  Arl{f(<);r}  =  i  J  (12) 

c— loo 

as  its  inversion  formula.  If  r“-1/(r)  is  absolutely  integrable  on  the  positive  real  axis  for 
some  a  >  0,  then  the  inversion  is  valid  for  c  >  a. 

The  Mellin  transform  of  the  product  of  r2  and  Poisson’s  equation  (1) 

+  r|: +  w)  Mr'*]  =  -7T5(r "  r')6(*  ~  (13) 

is  the  simple  form 

+  (14) 

This  fortunate  property  of  the  Mellin  transform  of  the  r  dependence  in  the  two-dimensional 
Laplacian  is  responsible  for  its  successful  application  [52]-[54]  to  potential  problems  in 
wedge-shaped  regions.  In  region  (2)  where  there  is  no  forcing  term,  this  procedure  gives 

(*’ +  w)  =  °-  (15) 

The  complex  variable  a  is  a  parameter  in  the  above  pair  of  ordinary  differential  equations 
in  4,  with  solutions 

¥a(s,  4>)  —  A(a)  sin(s^)  (0  <  4  <  a)  (16) 

{»,<!>)  =  B(s)sm(8<f>)  +  C(a)cos(s<f>)  (a  <  4  <  4')  (17) 

(a,  4)  =  D(a)  sin[s(^  -  jt)|  (4'  <  4  <  *)  (18) 


in  view  of  the  soft  boundary  conditions  (5)  and  (6).  Transformation  of  the  four  remaining 
conditions  (7)-(10)  gives  the  set  of  simultaneous  equations 


sin(sa) 

cos(sa) 

0 

0 


sin(sa)  cos(sa)  0 

—  cos(sa)  sin(sa)  0 

sin(s^')  cos(s^')  sin[«(7r  —  4')] 

cos (a4')  ~  sm(a4,)  -  cos [s(ir  —  4')} 


A(a) 

B(a) 


C(a) 

D(a) 


(19) 
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to  be  solved  for  the  coefficient  functions  in  (16)-(18).  Solution  of  these  yields  the  (soft) 
Mellin  transforms 


1  -  r  r,$ 

Ci  8 

■  sin[s(?r  —  0')]  sin(s^)/A(s) 

(0  <  <f>  <  a) 

(20) 

r" 

—  sinfs 

(ir  —  ^')]  {sin(s^)  +  T  sin[s(^  ■ 

-2a)]}/A(t) 

(a  <  </>  <  <f>') 

Ci« 

(21) 

*f  (»,*)- 

ru  .  r 
—  sin  [si 

(r  —  <f> )]  {sin(s^')  +  T  sin[a(^' 

-2a)]}/A(s) 

{<(>'  <  <t>  <  *) 

(22) 

with  denominator  function 


A(«)  =  sin(sTr)  +  T  sin[s(7r  —  2a)] 


and  dielectric  contrast  parameter 


€2  —  61 


*2  +  *1 

B.  Modifications  for  the  Case  of  Even  Symmetry. 

In  the  case  of  a  hard  ground  plane,  the  Neumann  boundary  conditions 


— *(r,  0)  =  0  (25) 

r)=0  (26) 

replace  the  Dirichlet  boundary  conditions  (5)  and  (6)  of  the  previous  section.  A  similar 

application  of  the  Mellin  transform  and  the  other  unchanged  boundary  conditions  yields 
the  (hard)  Mellin  transforms 

^(s,  <i> )  =  — — -  —  cos[s(tt  -  <f>')]  cos(s<t>)/A(s)  (0  <tj>  <  a)  (27) 

€x  S 
r'* 

(*>  =  —  cos[s(jt  -  <(>')]  {r  cos[s(^  -  2a)]  —  cos(s^)}  /A (a)  (a  <  <j>  <  </>') 

Cl®  (28) 

(*,  <t>)  =  — -  cos[s(»r  -  ^)]  {T cos[s(^'  -  2a)]  -  cos(s^')}  /A(s)  (<!>'<  4>  <  *0 

Cl*  (29) 

where  in  this  case  the  denominator  function  is 


A(«)  =  sin(«2r)  —  T  sin[s(ir  —  2a)].  (30) 

Note  that,  except  for  the  simple  scaling  by  6i  which  persists  from  the  original  source 
strength  chosen  in  (1),  the  presence  of  two  different  dielectrics  is  entirely  accounted  for  by 
T  in  all  of  the  above  transforms. 
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C.  Inverse  Mellin  Transforms  -  Preliminaries. 

The  zeros  of  the  denominator  functions  (23)  and  (30) 


A(s)  =  sin(sir)  ±  T  sin[«(jr  —  2a)] 


(31) 


are  central  to  the  Mellin  inversion  (12).  These  real,  simple  zeros  can  be  computed  via 
Muller’s  algorithm  [56]  for  arbitrary  half-angle  a,  or  solved  analytically  as  the  roots  of  a 
trigonometric  polynomial  when  a  is  a  rational  multiple  of  it.  This  procedure  is  demon¬ 
strated  for  the  particular  case  a  =  tt/3,  whereupon  the  variable  change  u  =  si r/3  in  (31) 
gives 

sin(3u)  ±  T  sin(u)  =  0  (32) 

which  factors  into 

[3  ±  T  -  4sin2(u)]  sin(u)  =  0.  (33) 

The  required  roots  are  now  explicitly  given  by 


sn 


+  3m 


(m  =  0,  ±1,  ±2, . . . ), 


(34) 


where  it  is  recalled  that  the  +r  (— r)  denotes  the  case  of  soft  (hard)  symmetry.  The  effect 
of  the  dielectric  material  (e2  ^  Ci  =>  T  ^  0)  on  the  potential  above  both  symmetry  planes 
is  a  regular  displacement  of  the  integer  poles  for  the  homogeneous  case  (c2  =  Cj  =►  T  =  0). 
The  index  n  in  (34)  is  a  denumerable  ordering  of  these  poles. 

D.  Inverse  Mellin  Transform  for  the  Case  of  Odd  Symmetry. 

As  r  — »  0  the  odd  potential  ip°  — » ►  0  and  the  dipole  behavior  ip°  ~  1/r  prevails  as  r  — ♦  oo 
in  the  far  field.  A  sufficient  Bromwich  contour  for  the  complex  integration  (12)  is  therefore 
guaranteed  for  the  choice  of  real  constant  0  <  c  <  1.  Complete  details  of  the  Mellin  inver¬ 
sion  for  the  potential  ^2(r,  <f>)  of  (4)  inside  the  dielectric  sector  are  provided,  whereupon 
the  final  forms  for  (r,  <f>)  and  tp*  (r,  <f>)  are  immediately  written  by  comparison. 

The  sin(s^)  factor  in  ^(s,  <f>)  of  (20)  together  with  the  transform  property  [52] 


A4_1  {sin(s^)F(s);  s}  =  [/(re'*)] 


(35) 


renders 

sin[s(7r  -  ^)](r7r)» 

{sin(ss-)  +  T  sin[s(7r  —  2a)]} 

G(.) 

the  desired  function.  The  integration  path  in  (36)  has  been  pushed  flush  against  the 
imaginary  axis,  and  the  principal  value  notation  invoked  to  properly  account  for  the  pole 
at  the  origin.  For  r  >  r',  closure  at  infinity  in  the  right-half  plane  gives  a  convenient 
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Fig.  3.  Bromwich  Contour  in  the  Complex  a— Plane 

contour  on  which  to  apply  Cauchy’s  integral  theorem,  and  is  shown  in  Fig.  3.  Let  Cn  be 
a  circle  of  vanishing  radius  p  — ♦  0,  centered  on  the  pole  an  with  residue 


lim  i  /  = 

p— *0  2tT«  Jcn  (*n) 


where  the  derivative  of  the  soft  denominator  function  is 

A'(s„)  =  Jrcos(s„ir)  4-  (jt  -  2o)r  cos[sn(7r  -  2a)].  (38) 

The  integral  around  the  origin 

lim  /  G(a)  da  =  ~  (39) 

p—q  2iri  Jc0  A'(0) 

follows  from  the  limit  sn  — »  0.  Let  positive  n  =  1,2, . . .  be  the  indices  of  the  poles  in  the 
right-half  plane,  and  let  negative  n  =  —  1,—  2, . ..  identify  the  poles  in  the  left-half  plane. 
The  odd  symmetry  of  the  function  A(s)  provides 

s_n  =  -a„  and  A'(-sn)  =  A'(sn).  (40) 

The  contribution  from  the  integral  around  the  infinite  semicircle  R  of  Fig.  3  is  zero  for 
r  >  r#,  while  closure  in  the  left-half  plane  is  appropriate  when  r  <  r'.  Cauchy’s  integral 
theorem  now  yields 

/w 5tl£fa(r/r')±‘’}  (r  s  r,)-  (41) 

which  together  with  (35)  gives 


4i(r, 4) 


_  r~  1  r  sin[sw(y -  ^)]8in(iw^)^ryt 


<i.A'(»n) 
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(r  %  r') 

for  the  sector  0  <  ^  <  a  (42) 


as  the  inverse  Mellin  transform  of  (20).  Since  all  three  of  the  transforms  (20)-(22)  are  of 
the  same  general  form,  the  remaining  two  held  expressions  are  apparently 


*)  =  =1  £  SSlhfe -  +1  (»in(»n^4-  rsin[,„(0  - Ml(r/f-)^..  (r  5 

n=l 

for  the  sector  a  <  <f>  <  <f>'  (43) 


and 


^,+  (r  ^  _  Zi  ”n^w^  ~  ^  {sin(*n^')  +  T ~  2q)]}  .yj 


>\±*r 


(r  g  /) 

for  the  sector  <f>'  <  <j>  <  7r.  (44) 


Note  the  reciprocity  in  the  last  two  expressions  for  <f>  «-►  <f>'.  The  above  reduce  to  the 
correct  potential  due  to  a  unit  line  charge  located  at  <j>’  =  ir/2  above  a  soft  ground  plane 
in  the  special  case  of  no  dielectric  wedge  (e?  =  a). 

E.  Inverse  Mellin  Transform  for  the  Case  of  Even  Symmetry. 

The  monopole  potential  lnr  due  to  the  original  line  source  and  the  in-phase  image 
is  the  dominant  feature  of  ipe  in  the  far  field  as  r  — ►  oo.  Therefore,  the  condition  of 
integrability  following  the  transform  pair  of  (11)  and  (12)  is  not  satisfied,  and  the  Mellin 
inversion  formula  cannot  be  directly  applied  to  the  functions  in  (27)-(29).  Mathematically, 
operating  with  the  ^-derivative  of  the  potential  temporarily  removes  this  troublesome 
logarithmic  variation,  which  is  then  restored  following  the  convergent  Mellin  inversion. 

The  ^—derivative  of  (27) 


d  «  / .  ^  1  ”  r  sin(s^)  u 

*)  -  — - ^ - r 


(45) 


is  transformed  in  the  manner  of  the  previous  section  to 

cos[*w(tt-  ^)]  sin(sw^) 


(r  $  r'l  (46) 


where  the  s„  are  now  the  hard  poles  in  (34)  and  — r  replaces  T  in  (38)  for  this  hard 
denominator  function  (31).  The  anti-derivative  with  respect  to  <f>  yields 


*M)  =  —  E 


e>  Si 


cos[*w(y  -  ^)]cos(<w^) 
snA'(s„) 


(r/r1)*-  +  -Mr)  (r  $  r') 

for  the  sector  0  <  <f>  <  a,  (47) 


where  the  boundedness  at  r  =  0  and  the  known  behavior  at  r  — ►  oc  specify 


{0,  for  r  <  r' 

— -ln(r/r'),  for  r  >  r' 

TTCi 


(48) 
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as  the  ^—independent  term  in  this  solution  of  Laplace’s  equation.  Similarly,  the  spatial 
potentials  from  the  related  transforms  (28)  and  (29)  are  thus 


rp-(r  =  “*[*„(*  -  ^)]  {r  cos[sB(^  -  2a)]  -  cos(s„^)} 

+ 1 v^(r)  (r  ^  r')  for  the  sector  o  <  <f>  <  <f>'  (49) 


and 


,  +  /  -1  y^COs[«n(lT-^)]{rcOs[sn(^ 

1  {)~~k  -.i'l'. 


-  2q)]  -  COs(Sn^)}  /_  ,_,s±tn 


) 


ir/r'Y 


+  tuT(r)  (r  ^  r')  for  the  sector  <j>'  <  <f>  <  ir.  (50) 


These  results  are  also  verified  for  the  special  case  of  a  unit  line  charge  located  at  <f>'  —  ir/2 
in  a  homogeneous  half-space  (ej  =  ei)  above  the  hard  symmetry  plane. 

F.  Results  for  the  Complete  Static  Solution. 

The  important  edge-behavior  of  the  dynamic  fields  in  the  immediate  vicinity  of  the  apex 
r  =  0  is  obtained  from  the  above  static  analysis.  A  quasi-static  philosophy  of  extracting 
the  dominant  behavior  of  a  solution  to  the  Helmholtz  (wave)  equation 

(V2  +  lb2)  1>(r )  =  0  (51) 

from  the  static  ( k  =  0)  solution  to  Laplace’s  equation 

V2V»(f )  =  0  (52) 

in  the  neighborhood  of  boundary  discontinuities  and  sources  is  well  known  and  successful 
([55],  [57],  and  [58]).  This  analytic  solution  for  the  static  wedge  problem  is  a  valuable 
resource  for  continuing  wave  studies. 

Contours  of  constant  potential  are  illustrated  in  Figs.  4-8,  all  for  the  case  of  wedge 
half-angle  a  =  tt/3.  Even  and  odd  potentials  are  computed  for  the  separate  boundary 
value  problems  and  combined  according  to  (3).  The  line  source  is  located  a  unit  distance 
(r'  as  1)  from  the  wedge  apex  at  the  (z,  y)  coordinate  origin.  The  domain  of  Figs.  4-7 
is  all  internal  to  the  unit  circle,  while  Fig.  8  encompasses  a  larger  area.  Unfortunately, 
the  poor  numerical  convergence  of  the  modal  series  at  r  as  r'  results  in  slight  kinks  in  the 
contours  near  this  unit  circle.  In  fact,  the  data  of  Fig.  8  is  already  smoothed  by  a  numerical 
averaging  based  on  the  mean  value  property  of  potential  solutions.  This  slight  blemish  in 
Fig.  8  has  no  effect  on  either  the  physical  interpretation  or  mathematical  integrity  of  this 
analysis.  Figs.  4-7  show  the  effect  of  varying  both  the  source  location  (<j>'  =  x/2  or  x)  and 
the  dielectric  constant  (e2  —  10  or  100)  of  the  wedge  relative  to  the  external  medium  of 
permittivity  ci  =  1.  When  £2  =  100  (Figs.  6  and  7),  the  dielectric  wedge  is  essentially 
an  equipotential  region,  which  is  expected  since  an  increasingly  dense  dielectric  (e2  — *■  oo) 
behaves  electrically  like  a  perfect  conductor. 
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Fig  .  5.  Contours  of  Constant  Potential  for  the  Dielectric  Wedge.  Case: 
r'  =  1,  <f>'  =  n,  a  =  ir/3,  =  1,  e2  =  10. 
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Fig.  6.  Contours  of  Constant  Potential  for  the  Dielectric  Wedge.  Case: 
r'  =  1,0'  =  ir/2,  a  =  x/3,  €1  =  1,  cj  =  100. 
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Fig  .  8.  Contours  of  Constant  Potential  for  the  Dielectric  Wedge.  Case: 
r'  =  !,<(>'  =  100  (deg),  a  =  ir/3,  €i  =  1,  c2  =  10. 


III.  Density  Contrast  Wedge  via  the  Kontorovich-Lebedev  Transform 
Following  D.  S.  Jones  [49],  the  Kontorovich-Lebedev  transform  pair 


OO 

9<y)  =  j  M«?\v)dy  (*►-«)  (1) 

0 

too 

*/(*)=  e !™0“2  J  eev\jv{x)g{t/)di/  (2) 

-too 

is  adopted  for  the  time-harmonic  analysis  (exp[+turt])  of  an  acoustic  plane  wave  incident 
upon  a  wedge  having  the  same  intrinsic  wave-speed  c  as  the  surrounding  medium,  but  with 
a  different  ambient  density  po.  Convergence  of  the  transform  is  assured  if 

l 

J I  f(v)  In  v\dy<oo  (3) 

o 


and  if 

OO 

J  f(v)v~^€'t  dy  <  oo  for  any  a  >  0.  (4) 

a 

The  geometry  is  still  that  of  Fig.  1,  but  the  acoustic  notation  is  adopted  in  this  section. 
Specifically,  a  line  source  of  unit  strength  is  located  at  the  cylindrical  coordinates  (r1,  <f>') 
in  the  external  medium  of  ambient  density  poi,  and  radiates  time-harmonic  waves  at 
the  radian  frequency  w.  A  wedge-shaped  scatterer  of  ambient  density  po2  occupies  the 
sector  —a  <  <f>  <  or.  The  key  to  applying  the  Kontorovich-Leoedev  transform  to  this 
idealized  transmission  problem  is  precisely  this  restriction  of  a  homogeneous  wave  speed 
c.  The  simplification  in  the  mathematics  is  obvious  in  view  of  the  ray  physics  at  planar 
interfaces  between  media  with  differing  density  but  with  the  same  wave  speed  and  therefore 
wavenumber  k  =  ufc.  This  observation  is  certainly  not  new,  as  evidenced  by  Chu’s 
[45]  elegant  derivation  of  the  impulse  response  of  this  density-contrast/isovelocity  wedge- 
scatterer. 

As  in  the  static  case,  the  mathematics  is  less  cumbersome  if  the  solution  is  split  into  its 
odd  and  even  components  with  respect  to  the  symmetry  plane  y  =  0.  In  particular,  the 
total  scalar  field  is 


rp{x,±y)=^e(x,y)±iff°(x,y)]  (y  >  0)  (5) 

where  the  odd  field  xp°  exists  above  the  soft  or  Dirichlet  ground  plane  and  the  even  field 
i>e  is  the  field  in  the  half-space  above  the  hard  or  Neumann  ground  plane  (see  also  Fig.  2). 


20 


A.  Odd  Symmetry  -  the  Soft  or  Dirichlet  Ground  Plane. 

A  statement  of  the  boundary  value  problem  consists  of  the  Helmholtz  equations 


(VJ  +  *!)  Mr*)  =  ~fa  -  r')H*  -  *') 

(a  <  <t>  <  ir) 

(6) 

and 

(vJ  +  *s)  Mr  *)  =  o 

(0<*< 

a) 

(7) 

together  with  the  Dirichlet  boundary  conditions  on 

the  soft  plane 

tMr)  0)  =  0  and 

V>i(r,7r)  = 

0, 

(8) 

and  continuity  of  both  the  pressure 

V»2(r,a)  =  V»i  (r, 

a) 

(9) 

and  the  normal  velocity 


— Tr^Hr>a) =  r- 

Pol  t)<P  P02  0<f> 


at  the  material  interface.  Additionally,  sufficient  decay  as  r  — ►  oo  (the  Sommerfeld  radia¬ 
tion  condition) 

lim  y/r  ( +  ikij>\  =  0  (11) 

r—oo  \0r  ) 

is  automatically  incorporated  into  the  chosen  integral  representation  by  the  Hankel  func¬ 
tion. 

Application  of  the  integral  operator 


OO 

J  drrH^\kr) 


to  the  partial  differential  equation  (6) 


r  Jr  (rJ^^r’^)  +JJ^2^i(r>^)-***Vi(r,^):B-p5(r-r'W--^')  (13) 

is  followed  by  two  consecutive  integrations  by  parts  to  remove  the  r— derivatives.  Explicitly, 
the  first  term  above  is  transformed  to 

OO 

J  Hl2)(kr)dr 

0 

oo 

=  »•  ~k  J  T^{rA)H^{kr)dr 


=0  by  (8)  and  (11) 


OO 

=  -*r^i(r,^)^2)'(*r)|  _Q+k  J  0i(r,^)  ^2)'(A:r)  +  krff^"(kr)^  dr. 


=0  by  (8)  and  (11) 
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Thus,  the  operation  (12)  on  (13)  gives 

OO 

/  V(2,*)  [(Jr)’ +  *rff<2>'(*r)  +  (tr)s«(2>(tr)] 

0  r  - * - - 

=i f2Hl2\kr)  by  BeMel’s  equation 

oo 

+  W  J  -  4'),  (15) 

0 

which  is  the  simple  differential  equation 

[w + "I) hM)  ~ (“>) 

in  terms  of  the  Kontorovich-Lebedev  transform 

OO 

AM)  =  JIL  ~  f  ^(tr)Mr^)dr.  (17) 

H*  \kr)i 

By  inspection,  the  source-free  equation  (7)  now  transforms  to 

+ ,,J) /*(‘'» = #  <18) 

with 

OO 

hM)  =  (2)1  /  \  ffl2)ikr)rj>7(r,<i>)dr .  (19) 

**»  \kr  )  j 

As  in  the  static  analysis  of  Section  II,  the  source  coordinate  <f>  =  <f>'  divides  the  external 
medium  (1)  into  two  source-free  subregions,  with  fields  (transforms)  distinguished  by  a  qF 


superscript: 

/a(i/ ,  <f>)  =  A(v)  sin(i /<j>)  (0  <  <f>  <  a)  (20) 

fi  4)  =  siB(v^)  +  C{u)  cos{vj)  (a  <  <f>  <  <t>')  (21) 

fti*,  =  D{v)  sin[i/(^  -  tt)]  <  <f>  <  it).  (22) 

Note  that  the  Dirichlet  boundary  conditions  (8)  have  been  used  in  writing  (20)  and  (22). 
In  the  source  neighborhood  <f>  =  the  function  f\  is  continuous 

St  O'* +')  =  /f  O'*  4>’)t  (23) 

leaving  the  derivative  operator  of  (16)  responsible  for  the  singularity,  which  upon  integra¬ 
tion  from  <j>'—  to  4>'+  yields 

/,+'(*',  - 1-  (24) 
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These  two  conditions  plus  the  transforms  of  (9)  and  (10)  result  in  four  simultaneous  equa¬ 
tions 


-sin(»/a)  sin(i/a)  cos(i/a)  0  "  r  q  ■ 

^•cos  (i/a)  —  cos(i/a)  sin(i/a)  0  B(u)  0 

P02  r,,  \  —  n 

0  sin(i/^')  cos(i/^')  sin[i/(jr  -  ^')]  ® 

0  cos(i — sin(t /<(>')  —  cos[i/(tt  —  ^')]  LX>(*/)J  ll/v. 


for  the  coefficient  functions  in  (20)-(22). 

The  above  is  necessarily  identical  in  form  to  the  matrix  equation  (19)  of  Section  II-A 
because  of  the  action  of  both  integral  transforms  on  a  single  partial  differential  operator 
in  both  media,  and  the  identical  boundary  conditions.  Therefore,  a  trivial  adjustment  in 
the  earlier  static  solution  produces  the  desired  transforms 

fi{y,  <t>)  =  ~  ---  sin[i/(ir  -  <j>')\  sin(i/^)/A(y)  (0  <  4  <  a)  (26) 

f\  4)  =  ~  sin[i/(jr  -  ^')]  {sin(i/^)  +  T  sin[i/(^  -  2a)]}  /A(i/)  (a  <  <j>  <  <j>') 


ft  4)  =  ~  sin[t/(ir  -  ^)]  {sin(x/<6')  +  T  sin[v(^'  -  2a)]}  /  A (i/)  {<j>'  <  <f>  <  r) 


with  denominator  function 


A(i/)  —  sin(i/x-)  +  T  sin[i/(ir  —  2a)] 


and  density  contrast  parameter 


P0\  —  P02 


PO\  +  P02 

The  Kontorovich- Lebedev  transform  pair  (1)  and  (2),  and  the  integral  definitions  of 
(17)  and  (19)  express  the  physical  fields  in  terms  of  complex  integrals,  such  as 

ioo 

V»2 (»*»  =  lim  f  etv,i/Jv{kr)Hl2){kr,)f2{i/,<l>)du  (0  <  <f>  <  a)  (31) 

—ioo 

and  similarly  for  V'f  (r,  4)  and  4>t  (r>  4)- 

B.  Even  Symmetry  -  the  Hard  or  Neumann  Ground  Plane. 

The  mathematical  analysis  of  the  hard  ground  plane  with 

•Orfair,  0)  =  0  and  (r,  ir)  =  0  (32) 


parallels  the  development  above  for  the  soft  case.  The  convergence  problem  (due  to  the 
logarithmic  static  behavior)  is  avoided  by  working  with  the  ^—derivative  of  the  wave 
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functions,  in  the  same  manner  as  in  Section  II-E.  In  this  case,  it  is  desirable  to  start  with 
the  derivative  of  the  field 

^(r,  ^)  =  (33) 

and  the  governing  partial  differential  equation  (in  region  (1),  for  example) 

~  (’•Ir&M))  =  -£«(•■  -  W  -  *')  (34) 

so  that  the  boundary  terms  that  result  from  the  integrations  by  parts  can  be  safely  set  to 
zero,  as  in  (14)  above.  This  procedure  is  carried  through  to  give  final  expressions  for  the 
transform  functions  of  (33)  which  must  be  evaluated  by  a  clever  combination  of  contour 
integration  and  possibly  some  asymptotics  to  extract  the  critical  physical  processes.  The 
static  analysis  of  Section  II  is  good  experience  and  background  to  continue  to  make  real 
progress  on  this  boundary  value  problem.  Enough  of  the  details  and  equations  are  pre¬ 
sented  here  to  enable  a  solid  continuation,  and  to  provide  a  basis  for  the  discussion  of  the 
impedance  boundary  condition  in  Section  IV. 

The  integral  operator  (12)  applied  to  (34)  yields 

(J^  +  -J)  SiK  *)  =  -«'(*  -  4f)  (35) 

in  terms  of  the  Kontorovich-Lebedev  transform 

OO 

*1 (",  *)  =  -;3 r—  /  *)  dr.  (36) 

Hv  \kr )  J 

The  subscript  and  superscript  notation  is  consistent  with  the  soft  case  above.  Physical 
boundary  conditions  on  the  ^—derivative  of  the  fields  can  be  directly  applied  to  the  trans¬ 
forms  g(u,  <f>),  but  boundary  conditions  on  ip  itself  involve  the  anti-derivative  of  the  g(v,  <f>) 
functions.  That  is,  the  Neumann  conditions  (32)  are 

52^,0)  =0  and  01- (»>,  *0  =  0,  (37) 


and  continuity  of  pressure  and  normal  velocity  at  the  media  interface  become 


a  a 


(38) 


and 


Poi  P 02 


(39) 


respectively.  Similarly,  the  boundary  conditions  due  to  the  impulsive  source  at  <f>  —  <f>'  are 


J  9i  (^  <P)d<f>  =  J  gi  (1 /,  <f>)  dtp  (40) 
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and 


(41) 


-  Pi  (*. ^')  =  — i* 

Evidently,  solutions  of  the  form 


<f>)  =  A(v)  sin(i/^)  (0  <  ^  <  a) 
g;  (j/,  4>)  =  B(v)  sin +  C(v)  cos(j/^)  (a  <  <f>  <  ^') 
0? (*>,  =  B{v)  sin[i/(<£  -  7r)]  (<t>'  <  <f>  <  rr) 


have  coefficient  functions  that  satisfy  the  linear  equations 


—  sin(t'a) 

P07 

cos(i/a) 

0 

0 


—  sin(va) 

—  COS^a) 
sin^^') 
cos(i/^') 


—  cos(ya) 

0 

‘ A{i ')' 

r°-| 

sin(i/a) 

0 

B{u) 

0 

cos(j/^') 

sin[i/(?r  —  ^')| 

C{v) 

1 

—  sin(t/^') 

—  cos[i/(?r  —  <j>')]  m 

lD(v)\ 

.0. 

(42) 

(43) 

(44) 


(45) 


Simple  formulas  for  these  transforms  are  also  available  through  comparison  with  the  ap¬ 
propriate  static  analysis,  which  in  this  hard  case  is  Section  II-E. 
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IV.  Impedance  Boundary  Condition 

A.  Homogeneous  or  Traditional  Impedance  Boundary. 

The  Leontovich  impedance  boundary  condition  avoids  any  consideration  of  the  held 
interior  to  the  wedge  region,  and  therefore  applies  in  the  case  of  a  highly  lossy  or  effectively 
impenetrable  wedge.  As  with  all  of  the  work  reported  here,  the  original  boundary  value 
problem  is  decomposed  into  odd  (soft)  and  even  (hard)  symmetry  components.  Only  the 
field  in  the  external  region  a  <  <f>  <ir  needs  to  be  examined  now,  so  the  subscript  notation 
is  abandoned,  but  the  superscript  is  kept  to  distinguish  the  fields  in  the  sectors  <f>^<(>'. 
The  impedance  boundary  condition  in  general  coordinates 

Q 

il>(r,a)  +  T)—il>(r,a ;)  =  0  (1) 

specifies  the  ratio  of  the  pressure  to  the  normal  velocity  at  the  material  boundary  <f>  — 
a.  The  numerical  value  of  this  surface  impedance  17  derives  from  the  actual  constitutive 
properties  of  the  dissipative  surface.  For  a  highly  lossy  material,  17  is  complex  with  a  phase 
angle  ir/4.  Note  that  the  limiting  values  17  =  0  and  17  — »  00  are  the  classical  Dirichlet 
and  Neumann  conditions,  respectively,  which  can  serve  as  partial  validation  cases.  The 
mathematical  treatment  of  the  even  symmetry  component  or  hard  ground  plane  problem 
(which  is  consistently  the  more  difficult  case)  is  illustrated  here. 

As  before,  a  line  source  of  unit  strength  is  located  at  (r',  <!>'),  in  the  space  between  the 
impedance  boundary  at  <f>  =z  a  and  ;he  hard  boundary  at  <f>  =  it.  The  forced  Helmholtz 
equation  is 

(Vs  +  *’)  V(r,  4)  =  -  i«(r  -  r')«  (*  -  *')  (.$*<.)  (2) 

where  the  two  physical  boundary  conditions  in  cylindrical  coordinates  are 

^(r,a)  +  ^^^(r,a)  =  0  (3) 

and 

~^(r,  tt)  =  0,  (4) 

plus  the  usual  Sommerfeld  radiation  condition. 

Let 

fri4>)=  (5) 

be  the  ^—derivative  of  the  field,  and  likewise  for  (2): 

li  «) +  hw({r'  *] + k'l(r'  = 46(r  -  -  «■  w 

Duplicating  the  progression  of  Section  III,  application  of  the  operator 

OO 

J  drrH^(kr)  (7) 

o 
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followed  by  integration  by  parts,  with  the  aid  of  the  boundary  conditions  and  Bessel’s 
equation  gives 

+  */2)  9^'  *)  =  ~6'^  ~  *')  W 


with  transform 


OO 


In  terms  of  ((r,  0),  the  impedance  condition  (3)  is 


a 

J  C  (»•,  <f>)  d<j>  +  (r,  a)  -  0. 


The  Kontorovich- Lebedev  transform  of  £(r,  </>)  is  denoted  alternately  by 
g~  ( v ,  <f>)  =  A(i/)  sin (i/^)  +  B{y)  cos(i/^)  (ot  <  <j>  <  <f>') 


=  C(i/)  sin[i/(^  -  *■)]  (^'  <  <t>  <  *),  (12) 

which  are  the  solutions  of  (8)  in  each  source-free  sector  on  either  side  of  the  source  discon¬ 
tinuity.  The  Bessel  recursion 

B?\kr)  _  B^^kr)  +  H^kr) 

kr  2u  [ 

and  the  anti-derivative  of  (11)  give  the  transformed  version  of  the  impedance  boundary 
condition 

~  {A[y  —  1)  sin[(i/  —  l)a]  +  B{y  —  1)  cos[(i/  —  l)a] 

+  A(i/  + 1)  sin[(i/  +  l)a]  +  B(v  -f  1)  cos[(v  +  l)or]} 

—  A(v)  cos(^o)  +  B(u)  sin(»/a)  =  0,  (14) 

which  is  a  linear  difference  equation  in  the  coefficient  functions  A( v)  and  B(v).  The  Dirac 
delta-function  at  <f>  =  <j>'  imparts  the  pair  of  boundary  conditions 

g~{v,  <!>')  = -\  (15) 


♦  9 

J g+{v,<l>)d<l>  =  J (16) 

or  equivalently 

C{y)  sin[i/(^'  —  7r)]  —  A(i/)  sin(i/^')  -  B{v)  cos(t/^')  =  —  1  (17) 

—  cos[i/(^<  —  *■)]  =  cosW)  +  ~~  sin(i/^').  (18) 

Although  apparently  this  is  not  the  same  approach  or  result  of  W.  E.  Williams  [7],  it 
is  interesting  and  significant  that  he  too  solves  a  difference  equation  in  modeling  the 
impenetrable  wedge. 
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B.  Inhomogeneous  or  Pseudo- Imp edance  Boundary. 

The  difference  equation  above  in  the  analysis  of  the  true  impedance  boundary  condition 
results  from  the  1/r  factor  in  (3).  In  the  perhaps  improbable  case  of  an  inhomogeneous 
surface  impedance  that  is  proportional  to  radius,  i.e.  rj  =  rrj1,  the  surface  condition  is 

^(r,a)  +  f7,^V’(r,a)  =  0  (19) 

and  no  shifting  occurs  in  the  argument  v  of  the  transform  coefficient  functions.  No  differ¬ 
ence  equation  need  be  considered,  and  the  set  of  linear  equations  (for  fixed  v) 

'  rfvsx&lyct)  —  cos(i/a)  if  v  cas(ya)  +  sin(i/a)  0  '  A(y)  0" 

cos(i/^')  -  sin[i/(^/  —  jt)]  B(  v)  —  1 

cos(i/^')  -  sin(i/^#)  -  cos[i/(^'  —  *)]  C(  v)  0 

(20) 

is  easily  solved. 

The  inverse  Kontorovich- Lebedev  transform  can  be  approached  from  an  analytical  ad¬ 
vantage,  and  an  inviting  opportunity  emerges  for  contour  deformation  and  asymptotics. 
Of  course  the  real  impedance  boundary  will  ultimately  yield  to  these  tactics,  but  in  future 
continuation  of  this  scattering  research.  For  the  present,  consider  these  far-field  (kr  — ►  oc) 
results 


tT(r,^)  = 


r  too 

«  - — itr  f  pj  ,kr>) {^cosW  -  a)l  -  sinW  -  a)]}  cos^'  -  tt)]  ^ 
2jrlfcr  J  *  y'u  sin[i/(7T  —  a)]  +  cos[i/(x  —  or)] 

—too 


(a  <  6  <  6,y\  (211 
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V.  Progress  Summary  and  Recommendations  for  Continuing  Research 


1.  The  complete,  analytical  solution  for  the  static  excitation  (Poisson’s  equation)  of 
the  dielectric  wedge  is  accomplished  using  the  Mellin  transform.  This  benchmark 
solution  is  believed  to  be  a  new  result,  and  is  a  valuable  tool  in  the  study  of  static 
edge  condition  behavior  at  the  apex  of  the  material  wedge. 

2.  The  Kontorovich- Lebedev  transform  has  been  successfully  applied  to  the  density 
contrast  wedge,  resulting  in  closed-form  expressions  for  the  integral  transform  rep¬ 
resentation  of  the  total  fields  in  both  the  interior  and  exterior  wedge  regions.  An 
asymptotic  evaluation  of  these  transforms  at  both  the  far  field  limit  (kr  — ►  oo)  and 
in  the  immediate  vicinity  of  the  apex  ( kr  -*  0)  is  of  great  interest  to  both  theo¬ 
retical  and  numerical  scattering  researchers.  Furthermore,  this  frequency  domain 
result  is  a  companion  to  Chu’s  [45]  transient  derivation. 

3.  The  difference  equation  present  in  the  transform  variable  of  the  Kontorovich- 
Lebedev  integral  needs  to  be  addressed  in  order  to  rigorously  solve  for  the  physical 
case  of  the  impedance  boundary  condition 


Such  a  rigorous  solution  provides  for  high  frequency  ray-launching  coefficients,  in 
the  manner  of  the  GTD  [10]  or  Wedge  Assemblage  [41].  Any  relationship  between 
these  promising  results  and  the  difference  equation  of  Williams  [7]  and  the  Mali- 
uzhinets  function  [6], [8-12]  must  be  affirmed  and  explained. 

4.  In  the  case  of  an  inhomogeneous  surface  impedance  that  is  proportional  to  radius, 
i.e.  r]  =  ri/,  the  above  Leontovich  boundary  condition  is 


*+4=°, 


resulting  in  closed-form  integral  transforms.  Inversion  to  spatial  coordinates  (r,  <f>) 
will  closely  follow  the  completed  Mellin  inversion  for  the  static  wedge.  If  this 
inhomogeneous  problem  is  deemed  to  be  of  physical  interest,  then  this  mathemat¬ 
ical  solution  is  of  great  value,  offering  great  physical  insight  into  the  asymptotic 
interactions  between  the  individual  scattering  mechanisms. 

5.  The  coupled  integral  equations  for  equivalent  surface  distributions  that  were  de¬ 
rived  during  the  initial  stages  of  the  present  research  could  form  the  basis  for  an 
alternate  method  of  attack.  However,  in  view  of  the  success  to  date  with  the 
Kontorovich-Lebedev  transform,  surface  integral  equations  will  not  likely  be  pur¬ 
sued  for  this  proposed  research. 

6.  Continuing  consultation  and  interaction  with  active  researchers  engaged  in  related 
electromagnetics  and  acoustics  work  should  be  maintained,  via  paper  submission, 
paper  reviewing,  telephone  and  written  correspondence,  the  1993  URSI  National 
Radio  Science  Meeting,  and  the  University  library  services. 
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Conclusions 


The  exact,  modal  series  solution  of  Laplace’s  equation  for  the  static  field  inside  and 
outside  of  the  dielectric  wedge  is  obtained  via  the  Mellin  transform,  and  is  readily  evalu¬ 
ated.  Similarly,  the  Kontorovich-Lebedev  transform  is  effectively  applied  to  the  problem 
of  time-harmonic  acoustic  interaction  with  the  density  contrast  wedge.  A  high-frequency 
asymptotic  analysis  of  this  transform  solution  is  a  promising  vehicle  with  which  to  un¬ 
derstand  and  extract  the  basic  physical  mechanisms  in  this  canonical  scattering  geometry. 
The  integral  transform  solution  of  the  wedge  having  impedance  boundaries  deserves  sub¬ 
stantially  more  mathematical  analysis. 
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